Fast detection of single-charge tunneling to a graphene quantum dot in a multi-level 

regime 
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In siiti-tunable radio-frequency charge detection is used for the determination of the tunneling 
rates into and out of a graphene single quantum dot connected to only one lead. An analytical model 
for calculating these rates in the multi-level tunneling regime is presented and found to correspond 
very well to our experimental observations. 
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Time-resolved charge detection on quantum dots is a 
powerful technique to straightforwardly extract single- 
particle transport properties such as the occupation prob- 
ability of quantum dots determined by the Fermi distri- 
bution function in the leads and the dot-lead tunneling 
rates Recording the charge-detector signal through 

radio-frequency (rf) reflection measurements Q can en- 
hance the time resolution drastically, enabling studies of 
systems with larger tunneling rates. 

To date, these types of experiments mainly study quan- 
tum dots in the single-level regime, and investigations 
involving transport through more than one level have fo- 
cused on the phenomenon of super-Poissonian noise 
Here, we present measurements and an analytical calcula- 
tion of multi- level tunneling rates using a graphene single 
quantum dot connected to one lead. 

Our charge-detection experiments were performed by 
incorporating a graphene nanoconstriction charge de- 
tector, capacitively coupled to a graphene single quan- 
tum dot, into an in site-tunable resonant circuit [7|. A 
schematic view of the latter and an atomic force mi- 
croscopy image of the former are presented in Figs. Ufa) 
and (b), respectively. The structure was fabricated 
through reactive ion etching on a mechanically exfoliated 
graphene flake [§[ . The parts of the graphene flake etched 
away are emphasized by dashed black lines in Fig. Qlb) 
to enhance visibility. The roughly 100 nm sized quantum 
dot (marked by the dashed red ellipse) and the nanocon- 
striction (designated by the single-sided arrow) can be 
tuned by in-plane gates and a global backgate. In spite 
of the symmetric design of the barriers forming the quan- 
tum dot, only the upper lead (source; represented by 
the double-sided arrow) was tunnel-coupled to the dot. 
Therefore, we were not able to measure direct transport 
through the graphene quantum dot. A more detailed 
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FIG. 1: (Color online) (a) Schematic representation of the 
circuitry used for rf and dc measurements. The dotted black 
capacitor is the total stray capacitance, and the dashed red 
boxes encase on-chip bias tees allowing for separation of low- 
and high-frequency signals. (b) Atomic force microscopy 
image of the graphene device studied in our experiments. 
Dashed black lines mark regions of graphene etched away to 
form a nanoconstriction at the location of the single-sided ar- 
row and a quantum dot encircled in dashed red. The nanocon- 
striction charge detector is attached to the resonant circuit as 
a resistive element R. (c) Conductance of the nanoconstric- 
tion measured via dc current (blue trace, left axis) and rf 
reflectometry (red trace, right axis). 



description of the fabrication procedure for this kind of 
devices can be found in Ref. |9| . 

Due to the large stray capacitance formed between the 
metallic bond pads and the highly doped silicon back- 
gate, high-quality rf matching of graphene samples us- 
ing the standard circuit design forming a resonance be- 
tween a series inductor and the total stray capacitance 
is difficult. One way to avoid this challenge is to re- 
move the backgate underneath the bond pads, which 
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FIG. 2: (Color online) (a) Time-dependent rf signal, 8-th 
order Bessel low-pass filtered at 200 kHz and sampled at 
500 kHz, revealing the times ti n and t out the quantum dot 
is occupied or unoccupied by a single negative excess charge, 
(b) Sketched energy diagram corresponding to the multi-level 
regime studied here. A quantum dot with a closely-spaced 
energy spectrum is connected to a lead with a tunnel rate F. 
The Fermi function of the lead is softened by a temperature 
T which is much larger than the level spacing AE but much 
smaller than the charging energy Ec- The difference between 
the Fermi level in the lead and the chemical potential of the 
quantum dot is denoted by A. 

however complicates the already involved fabrication of 
graphene nanostructures. We chose an alternative cir- 
cuit design consisting of a tunable capacitance in series 
to a parallel LCR circuit Q. Doing so and harnessing 
the large change in conductance of the nanoconstriction 
upon addition of a single charge to the quantum dot due 
to the proximity of the structures (AG ~ 0.08 x 2e 2 /h at 
G ~ 0.1 x 2e 2 /h) yields a very good charge sensitivity of 
5q « 3 X 10" 4 e/VSz 0,[TM1- 

All our measurements were performed in a variable- 
temperature insert at a temperature of roughly 2 K. 

The blue curve in Fig. QJc) shows the dc conductance 
(left axis) through the nanoconstriction charge sensor as 
a function of backgate voltage, while the red curve is the 
simultaneously acquired rf reflection coefficient. Both 
measurements reveal a transport gap in the backgate 
range between around —1 and 4 V. Subsequent measure- 
ments are performed in the regime of the dashed green 
box, in proximity of the charge-neutrality point of the 
graphene charge sensor [l4j. 

Tuning all gate voltages such that the quantum dot 
is at a charge-degeneracy point between N — 1 and 
N negative charges on the dot [H] and recording the 
charge-detector conductance via rf reflectometry allows 
for studying tunneling of single charge carriers in real- 
time, as shown in Fig. [2ja). This trace has been 8-th 
order Bessel low-pass filtered at 200 kHz - roughly two or- 
ders of magnitude faster than with previous time-resolved 
experiments on graphene quantum dots - and sam- 
pled at 500 kHz to ensure large enough signal-to-noise 
ratio for reliable electron counting From the aver- 
age time a negative charge spends inside (out of) the dot 
we can compute the dot occupation probability and the 
rate for tunneling out of (into) the quantum dot. 

At our elevated temperatures we have good reason [l8[ 
to assume that we are in a multi-level tunneling regime - 
an assumption that we will show to be correct in our dis- 



cussion of the experimental data. Thus we now calculate 
the tunneling rates of a quantum dot connected to one 
lead in the multi-level regime schematized in Fig. HTb). 
Our analysis closely follows the work of Beenakker [l9| 
in which all the ingredients for the subsequent analysis 
are provided, and we also adopt the notation thereof. 

The multi-level regime of Coulomb blockade is char- 
acterized by a temperature k^T much larger than the 
single-level spacing AE but much smaller than the charg- 
ing energy E c 

A£«feT«% (1) 

The latter inequality is validated by the fact that we ex- 
tracted an approximate charging energy of 10 — 20 meV 
typical for graphene quantum dots of this size (see for 
instance Ref. [3), whereas we assume the former to be 
true since we have not observed any sign of excited states 
in either the similarly sized charge detector or the dot it- 
self. Nevertheless, we shall assert the inequality's validity 
later-on. 

It has been shown in Ref. [19] that the rate for tun- 
neling into a quantum dot occupied by N — 1 charges is 
given by the sum over all quantum dot states p of their 
respective couplings T p multiplied by the probability to 
have a state at energy E p available in the lead (given by 
the Fermi distribution function /) and the conditional 
probability 1 — F eq {E p \N — 1) of having state p empty in 
an N — 1-fold occupied dot in equilibrium 

Tin =E P T p [l~F cq {E p \N-l)} 

xf{E p + U N -U N -i-E F ), (2) 

where Un is the electrostatic energy of a dot contain- 
ing N negative charge carriers, and Ep is the chemical 
potential in the lead(s). Beenakker has also elucidated 
that F oq {E p \N) can be expressed as a Fermi function 
f(E p — /in) in the high-temperature limit, where the 
chemical potential ^n of the dot is determined by 

oo 

J2f( E P~^)= N - ( 3 ) 

p=0 

By assuming the couplings T p to be equal and inde- 
pendent of energy, we can replace the sum by an integral 
over energy times density of states p which we will also 
assume to be energy- independent. Thus we can perform 
the integration which yields 

_ pT [~Un + Un-i + E F - hn-i] 

1 - exp [{Un - Un-i - E F + HN-i)/kB,T\ 



~ H 1 - exp(-A/fc B T) ' v ; 

where A is the energy detuning between the chemical po- 
tentials in the lead and the dot (note that our definition 
of A has an additional minus sign compared to the one 
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of Ref . [19[ ) . Analogously, we can compute 

Tout =E P T p F cq (E p \N) 

x [1 - f(E p + U N - Efo-i - E F )] 

= PT l-exp(A/fe B T)- (5) 

In this last equality we have made use of the fact that at 
our relevant energies /j>n-i ~ /Mr- 

Adding these two expressions leads to the inverse of 
the correlation time of the random telegraph signal [13] 
given by 

7 = r in + r out = rpA — — -. (6) 

cosh(A/KBJ J — 1 

We want to emphasize that in contrast to the above ex- 
pression, the sum of tunnel rates is constant with respect 
to detuning in the single-level regime if the couplings for 
tunncling-in and -out are equal (and are only varying 
over a range of the order of k^T otherwise). Therefore, a 
strongly varying sum of tunneling rates excludes a single- 
level regime. 

The expected rate of events te (the number of times 
per second an electron tunnels in and out of the quantum 
dot) can now be calculated to be 



Tin -To 
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2 sinh( A/ k B T) 

— fc B Tcoshr 2 
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(7) 



equivalent to what Beenakker found for the conductance 
through a dot in the multi-level regime (and what Kulik 
and Shekhter found for transport through granular met- 
als plj ) since in all instances, tunneling-in and -out are 
treated on an equal footing. 

Equipped with these analytical expressions for the tun- 
neling rates, we turn to our experimental results. Figure 
[3](a) shows the fraction of time the quantum dot is de- 
void of the AT-th excess charge when sweeping over the 
{N — 1, N} charge-degeneracy point, reflecting the Fermi 
distribution function of the lead. Fitting the measured 
data (blue dots) with a Fermi function (solid black curve) 
results in a temperature of T — 1.97 ± 0.03 K. Using 
this temperature, we can compare the number of tun- 
neling events extracted from the same time traces with 
the expected number in a single- (dashed green line; see 
Ref. [19|]) and multi-level (solid red line; Eq. [7]) regime 
in Fig. EJb). We find that the solid red line correspond- 
ing to multi-level tunneling matches the measurements 
better [22j]. 

In Fig. EJc) we show the experimentally extracted 
rates for tunneling-in (blue circles) and tunneling-out 
(red crosses). Additionally, we plot the analytical func- 
tions for the rates given by Eqs. 0] and [3] as solid lines. 
The only fit parameter for these curves is the product 
Tp of tunnel coupling and density of states which was 
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FIG. 3: (Color online) (a) Electronic occupancy of the lead 
as a function of dot-lead energy detuning A. Every blue data 
point is the fraction of time the dot is occupied by TV — 1 
(instead of N) negative charges as extracted from traces as 
shown in Fig. [2ja). The solid black line is a fit to a Fermi 
distribution function in the lead, yielding a temperature of 
T = 1.97 ± 0.03 K. The dashed black lines indicate the mag- 
nitude of temperature, (b) Number of tunneling events in a 
10 ms long trace. The dashed green and solid red lines are cal- 
culated rates for a single and multi- level regime, respectively, 
using the temperature extracted in (a) and the height of the 
data points, (c) Tunneling rates determined through the av- 
erage dwell times in (red crosses) and out (blue circles) of the 
quantum dot. The solid lines are calculations using Eqs.|4]and 
[5] in the main text. For sake of clarity, the statistical error 
bars are only shown for one in 20 traces, (d) Measured (blue 
squares) and analytically determined (Eq. [6] solid red line) 
sum of the multi-level tunnel rates shown in (c). The dotted 
green line indicates the behavior in the single-level regime. 



chosen such that the experimental results coincide with 
the calculations at zero detuning. The error bars indicat- 
ing the uncertainty of the determination of the tunneling 
rates from the dwell times [|[ are only shown for one in 
20 data points to enhance visibility. Finally, the sum of 
these tunneling rates is presented in Fig. [3Jd) for both, 
experiment (blue squares) and single- (dotted green line) 
as well as multi-level (solid red line; Eq. [6]) models. As 
we observe a linear increase in the sum of the rates for 
detunings above k^T, we can be certain to be in a multi- 
level regime. 

In all plots of Fig. [31 the horizontal axis has been con- 
verted to energy using a lever arm set via the relative in- 
fluence of the swept backgate and the dot source lead [23[ . 

The qualitative and quantitative agreement between 
our experiments and theoretical calculations is striking, 
especially considering that only one free parameter Tp is 
available, since temperature has been fixed by the Fermi 
fit in (a). Nevertheless, we observe a deviation from the 
model for large detuning, where the experimentally ob- 
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tained values lie distinctly above our expectations. This 
may originate from the fact that, in contrast to our as- 
sumptions, either the density of states p or the tunnel 
couplings T p in the dot are not entirely independent of 
energy. Furthermore, the couplings may be different for 
each state p. Also, it is conceivable that only a finite num- 
ber of states contribute to tunneling. From the energy- 
dependence of the sum of tunneling rates in Fig. [U[d), 
though, we can definitely exclude single-level tunneling. 
Missed events due to our large yet finite bandwidth lead 
to a further enhancement of the numeric value of the 
measured tunnel rates for large A and can therefore not 
explain the deviations. 

In conclusion we have performed time-resolved charge 
detection measurements on a graphene quantum dot us- 



ing rf reflectometry with large bandwidth. From these 
experiments we have extracted the rates for tunneling-in 
and tunneling-out of the dot, which are in good agree- 
ment with our analytical model for multi-level tunneling 
to a single lead. This type of fast charge detection on 
graphene quantum dots opens possibilities for studying 
interactions of Dirac Fermions in the framework of full 
counting statistics and is a great instrument for deter- 
mining spin relaxation times in graphene quantum dot 
circuits. 
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